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An experimental investigation of a low power arcjet plume was conducted using emission spectroscopy. A 
laboratory -model arcjet incorporating a segmented-anode was run on simulated hydrazine at a flow rate of 5 x10 s 
kgls. The complete visible spectrum measured in the exit plane of the arcjet showed the presence of N 2 , N 2 +, NH, 
and H. Radial intensity profiles for the H a , H p, and the NH A 3 11 X 3 L (0,0) transitions at four different axial 
locations were measured. These line of sight intensity measurements , spaced 0.05 mm apart , were deconvoluted to 
give the radial intensity distribution using an inverse Abel transformation. The ratio between the intensities from 
the H a and Hp transitions indicated a non-Boltzmann energy distribution between excited states in the plume. Axial 
intensity profiles taken on center line indicated the decay rate of excited states in the plume . An electron number 
density of 2xl0 13 cm' 3 at the exit plane was determined based on Stark broadening of the Hp line. Rotational 
temperatures of 750 K, 1 750 K, and 2500 K were determined for N 2 , N 2 +, and NH respectively. The results of this 

investigation demonstrate that the location of the current attachment on the anode has a measurable effect on the 
electronically excited species in the plume and that dissociation is the dominant frozen flow loss mechanism in low 
power arcjets. 


Introduction 

Research and development on one kilowatt-class 
arcjets to be used for north-south station keeping of 
geosynchronous satellites has been underway for 
several years. 110 As a result, this propulsion system 
has rapidly evolved from an unsteady laboratory 
device to a flight type unit. An extended duration 
lifetest has shown thruster lifetimes greater than 1000 
hrs. 2 The power processing unit and other necessary 
components have also been developed to a flight ready 
status. 3 Recent research efforts have been focused on 
addressing the remaining user concerns, answering 
fundamental questions about arcjet operation, and 
improving and extending the capabilities of the arcjet. 
One primary concern for arcjet system integration 
onto a communication satellite is the exhaust plume. 
Unlike the resistojet, which also heats decomposed 
hydrazine and expands it through a converging- 
diverging nozzle, the arcjet heats the propellant to a 
greater extent in an arc discharge. The resulting 
plume is therefore more energetic than a resistojet 
plume. As a result of the arc, electronic excitation, 


ionization, and dissociation occur in the nozzle. 
Immediately downstream of the nozzle, the plume is a 
rapidly-expanding, partially-ionized plasma in which 
recombination and de-excitation take place. In order 
to quantify the effects of the plume on different 
spacecraft subsystems, various experimental and 
analytical capabilities have been employed. 
Experimental investigations, such as those conducted 
by Carney 4 ’ 5 using electrostatic probes, have provided 
data on the electron number density and the electron 
temperature in the far field of the plume. The effect 
of the arcjet plume on the transmission and reception 
of communication signals has been modeled based on 
these measurements. 6 Litde, however, has been done 
to characterize the near field plume. 

Many studies have been conducted to investigate 
arcjet operating characteristics. Experiments have 
demonstrated the relative performance of different 
geometries, 7 electrical configurations, 8 and operating 
conditions. 9 These empirical studies have given 


insight into the physics governing the operation of 
the arcjet. However, information such as arc 
temperature and the nature of the arc attachment on 
the anode are not well known. Such data would 
permit a more thorough understanding of the internal 
energy deposition processes which dictate the 
operation of the arcjet. Information about the near 
field plume will provide information about 
phenomena occurring within the arcjet and result in a 
more thorough understanding of the plume which is 
needed for an accurate impacts assessment 

This report details an experimental program 
initiated to characterize the near field of the arcjet 
plume. Measurements were taken in the near field of 
a laboratory-model arcjet with a segmented anode 
using emission spectroscopy. The effect of anode 
current distribution on the plume was investigated. 
These results, in conjunction with the results from a 
companion paper, 10 provided information about the 
energy loss mechanisms in kilowatt-class arcjets. 


Nomenclature 

A transition probability 

a curve fit parameter 

b curve fit parameter 

B € rotational constant 

B v rotational constant for vibrational state 
c velocity of light, 3.0xl0 8 m/s 

C tm constant, (eqn 7) 

E h ionization energy for hydrogen, 

/ variable in transformed domain 

G Fourier transform of I 

& discretized values of G 

g statistical weight 

h Planck’s constant, 6.626x1 O' 34 J s 

/ intensity, W 

/ discretized values of I 

1 index for summation 

J rotational quantum number 

; VT 

J 0 zero order Bessel function 

k Boltzmann's constant, 1.380x1 O’ 23 J/K 

/ index for summation 

m species mass, g 

N number of discrete points 

N € electron number density, cm 3 

n quantum number 

Q r partition function, (eqn 8) 

r radial distance, cm 

T temperature 

T gx excitation temperature, K 

T e electron temperature, K 

T rol rotational temperature, K 

v vibrational quantum number 

2 species charge + 1 

a t rotational constant 


cl n 2 theoretical FWHM (Ref. 12) 

e emission coefficient 

AX D h 2 Doppler FWHM, A 

AX SIJ2 Stark FWHM, A 

A wavelength, A 

A 0 wavelength of line origin, A 

\x species mass, AMU 

v wave number, cm* 1 

v 0 wave number of molecular band head, cm* 1 

superscripts: 

" upper excited state 

lower excited state 
subscripts; 

1 upper excited state 

2 lower excited state 

Analysis 

Excitation temperature. Within a plasma, 
energetic collisions can result in electronic excitation. 
The population of an electronically excited state is 
proportional to the population of lower energy states, 
the electron number density, the electron energy 
distribution, and the appropriate collisional cross 
sections. As the electronically excited species relax, 
quantized energy is spontaneously emitted. The 
emitted energy from this radiative de-excitation is 
proportional to the population of the excited state and 
the probability that the electronic transition will 
occur. In an optically thin plasma, the emitted 
energy passes through the plasma without being 
reabsorbed. When the energy distribution between 
electronically excited states can be described by a 
Boltzmann equilibrium distribution, and the plasma is 
optically thin, the ratio between the energy emitted 
from two different transitions can be written as: 

li = exp(- (E r E 2 )/kT ot ) (1) 

2 A 2g2 X, 

From this expression the excitation temperature, 
T cx , can be calculated. The significance of this 
temperature is dependent on the degree of equilibrium 
within the plasma. For an optically thin plasma, if 
collisional processes occur much more rapidly than 
radiative processes, the plasma is considered to be in 
local thermodynamic equilibrium (LTE). For a 
plasma in LTE, any point can be described by the 
local values for temperature, density, and chemical 
composition. However, in many cases the conditions 
for LTE are not met. In a non-LTE plasma all 
equilibrating collisional processes are not 
significantly faster than other plasma processes. As a 
result, the energy distribution for different species can 
not be described by a single temperature and multiple 
temperatures are employed. Atoms may have an 
excitation temperature and a kinetic temperature. 
Molecules may have several different temperatures 
including the rotational temperature, vibrational 
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temperature, excitation temperature, and kinetic 
temperature. In this study, the H a (6562.8 A) and 
Hp(4861.3 A) electronic transitions from the Balmer 
series were used to calculate an excitation temperature 
based on the line ratio method presented above. 


The suitability of the assumption of a Boltzmann 
energy distribution between the third and fourth 
energy levels of hydrogen, ie. the upper excited states 
of Hq and Hp, can be estimated using the criteria for 
partial local thermodynamic equilibrium (PLTE) 
presented by Griem, see Ref. 11. For an optically 
thin plasma, this requires the electronic collisional 
rates for the upper excited states to be at least ten 
times faster than the radiative decay rates. The 
required electron number density for equilibrium 
between electrons and upper excited states in a 
hydrogenic plasma depends on the electron 
temperature. 


N c > 7x10 


17/2 


kT 
z E h 


1/2 


( 2 ) 


Electron number density. Spectral line 
emissions are broadened by several different 
mechanisms. These include Doppler broadening and 
Stark broadening. Instrumental limitations can also 
lead to apparent broadening. Doppler broadening is 
the result of a frequency shift of the emitted light due 
to the motion of the radiating panicle in relation to 
the observer. Because the radiating species have some 
relative statistical velocity distribution, not all the 
emitted light is shifted the same amount and the 
spectral line is broadened. This relative velocity 
distribution describes the kinetic temperature of the 
emitting species. For a Maxwellian velocity 
distribution, the line shape due to Doppler broadening 
will be Gaussian and the full width at half maximum 
(FWHM) for the line can be written as: 

2 tfl ... 

A X di /2 = [(8kT In 2)/(mc ) ] A.q (3) 

Doppler line widths for several different temperatures 
are shown in Table I for the Hp line. 


Table I: Calculated FWHM due to Doppler 
Broadening 


T.K 

A^Dl/2’ A 

1000 

0.078 

2000 

0.110 

3000 

0.135 

4000 

0.156 

5000 

0.174 

6000 

0.191 

7000 

0.206 

8000 

0.220 

9000 

0.233 

10000 

0.246 


Stark broadening produces Lorentzian line shapes 
and is caused by the interactions of the emitters with 
charged particles. Theories have been developed to 
predict how ions or electrons affect a line profile. 12 
The functional dependency of the electron number 
density on the FWHM can be simplified: 12 

-9 2/5 

A A si/ 2 = 2.50x10 ctj/^Ng (4) 

This expression was used to predict the FWHM for 
the Hp line as a function of electron number density. 
Calculated values are displayed in Table II. 

Table II: Calculated FWHM due to Stark 


Broadening 


N e , cm' 1 

A^Sl/2> A 

lxlO 12 

5.56x1 O’ 3 

5xl0 12 

4.58xl0' 2 

lxlO 13 

1.14x10 1 

5xl0 13 

9.38X10 1 

lxlO 14 

2.33x10° 


The combined effects of Doppler and Stark 
broadening can be described using a Voigt profile. 12 
If the kinetic temperature is known, the contribution 
to the Voigt profile due to Doppler broadening can be 
estimated. From the ratio between the observed line 
width and the estimated Doppler line width, the 
amount the line was Stark broadened can be 
determined. This gives the electron number density. 

Rotational Temperature. The different 
vibrational and rotational molecular energy states give 
rise to vibrational-rotational band spectra. During the 
transition from a higher electronic energy level to a 
lower electronic energy level, the vibrational energy 
state may change. Additionally, within any 
vibrational state there are many discrete rotational 
states. These different rotational states, separated by 
small energy differences, result in the many individual 
lines characteristic of the molecular spectra. From a 
rotating vibrator approximation, the wave numbers 
for the different rotational lines for the three different 
branches can be estimated via the following 
equations: 13 

v =v 0 + 2B , v +(3BVB"v)r+(B , vB\)r 2 (5a) 

v =v 0 + (B' v -B" v ) J , ’+(B , V -B" v ) J" 2 (5b) 

V = v o - (B 1 V +B" v ) J"+(B' V -B" v ) J" 2 (5c) 

In these, B v is a vibrational constant that varies with 
the vibrational quantum number. 

B v = B e -<x e (v+i) (6) 

Assuming a constant transition probability for 
all the lines of a rotational-vibrational band and a 
Boltzmann energy distribution between the different 
rotational levels within a given vibrational stale, the 
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intensity of each of the rotational lines can be 
calculated, 

4 

I = (J'+J"+l) expjBV (J'+ 1) hcAT] (7) 

Qr 

where C em is a constant and the partition function Q 
can be approximated by 


Using equations (5) - (8) the complete molecular 
spectra for a particular transition can be calculated. 
The temperature used to determine the intensity 
defines the distribution of rotational energies in the 
vibrational state and will be referred to as the 
rotational temperature. The fact that the rotational 
temperature determines the intensity distribution in a 
rotational-vibrational band was used to determine the 
rotational temperature from the measured spectra. For 
this an iterative process was used to determine the 
temperature giving the best agreement between the 
calculated band and the measured band. 

Apparatus 

An experimental system was designed and built 
to take spectroscopic measurements in the plume of 
kilowau-class arcjets. The arcjet was installed within 
a cylindrical vacuum chamber measuring 1.5 m in 
diameter by 5 meters in length. The chamber's 
pumping system consisted of four 0.82 m oil 
diffusion pumps (each rated at 30,000 1/s), a lobe type 
blower, and two piston type roughing pumps. The 
ambient pressure during arcjet operation was 
nominally 0.067 Pa. 

The arcjet Fired along the major axis of the tank 
and was mounted on a two-axis, computer-controlled, 
linear actuator system. This actuator system was 
capable of moving the thruster vertically and 
horizontally relative to the collection optics with a 
positional accuracy of about 0.025 mm. Both 
actuators were mounted on a single platform 
vibrationally isolated from the vacuum chamber. 

The optical system, shown in Figure 1, was 
designed to image the side view of the plume on the 
entrance slit of a monochrometer. The overall 
magnification of the system was 2.8. The primary 
collector was a 150 mm diameter parabolic mirror 
with a 406 mm focal length. This mirror was 
positioned 1,4 m from the center line of the plume, 
resulting in a solid angle of 8.9xl0 3 steradians. A 
second parabolic mirror with a diameter of 100 mm 
and a focal length of 228 mm was used to provide the 
desired magnification and focus the image of the 
plume on the entrance slit of a 0.64 m Czemy-Tumer 
type monochrometer. Data were taken with the 
entrance slit of the spectrometer opened to a height of 
3 mm and a width of 10 microns. Consequently, the 
cross sectional dimensions of the volume of the 


plume measured at each point was approximately 1 
mm in height by 3.6 microns in width. The 
spectrometer had a 110 mm x 110 mm holographic 
ruled grating with 2400 grooves/mm. This grating 
was blazed for use between 2500 A and 7500A. The 
theoretical linear dispersion of this spectrometer 
system was 0.06 A/mm. 

Two different photodetectors were used. A side-on 
photo-multiplier tube (PMT) with a quartz window 
and a multi-alkali photocathode was mounted on the 
axial exit port of the spectrometer, and a 1025 
element intensified linear photodiode array was 
mounted on the lateral exit port. Both detectors were 
therm oelectrically cooled to minimize the dark current 



Figure 1. Experimental Apparatus and Optical 
Configuration, Not to Scale 

The arcjet used in this investigation was a 
laboratory-model thruster with a segmented 
anode/nozzle. The thruster, which has been used in 
other investigations, has been described in detail 
elsewhere. 10 Briefly, the major components include a 
2% thoriated tungsten cathode 3.2 mm in diameter 
terminated in a cone with a 30 degree half angle 
(which corresponds to the convergence half angle of 
the anode). The propellant was injected tangentially 
through an annular injection disk to aid in arc 
stabilization. 

A radiation-cooled, segmented anode/nozzle 
consisting of alternating conducting and 
nonconducting segments was used. A simplified 
cross-sectional diagram of both electrodes is shown in 
Figure 2. The size and contour of the anode were 
chosen to approximate a standard laboratory-model 
configuration. The nozzle diverged with a 20 degree 
half angle, had constrictor dimensions of 0.64 mm in 
diameter by 0.25 mm in length, and an area ratio of 
225:1. The constrictor segments were numbered one 
through five. Segment number one corresponded to 
the furthermost upstream segment which included the 
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constrictor. Segment number five corresponded to the 
last downstream segment. Segments one and five 
were fabricated from 2% thoriated tungsten, 
Segments two t three, and four were fabricated from 
tantalum. The nonconducting insulators between the 
segments were fabricated from boron nitride. The 
entire anode was held in compression with four bolts, 
each with a spring to accommodate the thermal 
expansion of the anode. 


^-rmr — , 


Segment 1 2 3 4 5 


[^Tungsten 
nTantulum 
| Boron Nitride 


Figure 2: Simplified Schematic of Electrode 
Configuration 


linear diode array. This permitted calculation of 
accurate intensity ratios for the determination of the 
excitation temperatures from H a and Hp. A tungsten 
ribbon lamp was used as the calibration light source. 
The lamp was placed inside the chamber at the same 
location as the arcjet. The temperature of the 
filament was measured with a two color pyrometer. 
Simultaneously, the intensity of the incident radiation 
was measured at the wavelengths of interest using the 
diode array. The energy flux through the entrance 
slit of the spectrometer at a particular wavelength 
was then calculated using Planck's spectral 
distribution of emissive power for gray bodies. The 
emissivity used in this calculation was from a 
temperature and wavelength dependent polynomial 
expression from Larrabee. 14 The comparison of the 
measured values with the calculated energy flux 
resulted in a calibration of the measured intensity 
versus actual intensity. 


The electrical configuration of the anode was 
such that any of the segments could be manually 
switched in or out of the circuit. The current through 
each segment was measured with a Hall effect current 
probe. The voltage potential between the cathode and 
each segment was measured using a digital 
multimeter. The arcjet was run with a pulse-width- 
modulated, laboratory-model power supply. Arc 
breakdown was initiated by a 3 kV, microsecond 
duration starting pulse. 

All tests were conducted using a propellant 
mixture of hydrogen and nitrogen in a 2:1 ratio. This 
propellant mixture simulated the decomposition 
products of hydrazine. The flow rate of each gas was 
regulated with a mass flow controller. 

Procedure 

For all tests the arcjet was run with a total 
propellant flow rale of 5.0 x 10* 5 kg/s at a total 
current of 10 Amperes. These are typical conditions 
anticipated for flight application. Prior to each series 
of experiments the flow meters and the current probes 
were calibrated. In addition, the power supply was 
operated with a resistive load to verify the 10 Ampere 
output. Three different electrical configurations were 
chosen to minimize the total number of experimental 
runs due to the limited durability of the individual 
anode segments. The electrical configuration for each 
of these three configurations is shown in Table III. 


Table HI: Electrical Configurations 


Segment 

i 

2 

3 

4 

5 

Configuration 1 

0 

X 

X 

X 

X 

Configuration 2 

O 

0 

O 

o 

X 

Configuration 3 

X 

X 

X 

X 

o 

X denotes conducting, O denotes nonconducting 


After the intensity calibration was completed, the 
thruster was installed on the actuator system and 
realigned relative to the optics. A small incandescent 
light bulb was placed inside the anode such that a 
portion of the filament was visible from the side 
view. The thruster was then positioned horizontally 
and vertically relative to the primary collector until 
the image of the light bulb filament was properly 
aligned on the entrance slit. At this time the second 
parabolic mirror was adjusted such that the image was 
properly focused. The light was then removed and the 
chamber was pumped down to the high vacuum 
condition. 

Prior to taking spectroscopic data the arcjet was 
run until steady-state conditions were obtained. This 
was determined based on constant voltage/current 
characteristics, constant anode/nozzle temperature as 
measured with a two-color pyrometer, and negligible 
additional movement of the nozzle exit plane due to 
thermal expansion. Steady-state anode temperatures of 
approximately 1460 K were measured on the outer 
surface. Thermal expansion of the arcjet was 
monitored optically using the spectroscopic system. 
The arcjet was positioned such that the nozzle exit 
plane coincided with the entrance slit of the 
spectrometer. However, as the temperature of the 
arcjet approached steady state, the arejet's length 
increased due to thermal expansion and the nozzle exit 
plane was no longer imaged on the entrance slit. Due 
to limited conductive paths between the segmented 
anode and the anode housing it took approximately 30 
minutes for this equilibrium condition to be satisfied. 
The thruster was then repositioned axially until the 
nozzle exit plane was again properly imaged on the 
spectrometer’s entrance slit. 


Before pump down of the vacuum chamber an 
absolute intensity calibration was performed using the 


Intensity measurements were then taken with 
either the diode array or the photomultiplier tube. The 
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intensity measurements taken with the diode array 
consisted of radial and axial scans through the plume 
for a given spectral region. Radial scans were taken 
at four different axial locations: the exit plane, and 
2.5, 5.0, and 7.5 mm downstream of the exit plane. 
Each radial scan consisted of 257 spectra separated by 
0.05 mm. The axial scans were taken along the 
plume centerline. Each axial scan consisted of 200 
spectra, each separated by 0.05 mm. For each of the 
axial and radial scans the spectral regions containing 
the H a and Hp emissions were monitored. Radial and 
axial scans were also taken for NH A 3 n -> X 3 I (0,0) 
and (1,1) transitions. One radial scan was taken at the 
nozzle exit plane for N 2 + B 2 !^ X 2 I g (0,0) and (1 ,1) 
transitions. In most cases total intensities were 
measured by integrating over an interval two hundred 
times the full width at half maximum of any 
particular line as discussed in Ref 12. The dark 
current signal from the detector was subtracted from 
the measured total intensity as was the signal due to 
continuum radiation. The uncertainty in the 
intensities measured with the linear diode array were 
found to be within ±5 percent based on observed 
temporal fluctuations. The data taken with the 
photomultiplier tube consisted of wavelength scans at 
one axial location on the plume centerline. These 
data covered the region from 3200A to 7200A. Data 
were taken at 0.25 A increments and the intensity 
signal was integrated over 0.5 seconds. 

As previously indicated, both the axial and radial 
scan data were taken using the linear diode array. A 
computer controlled data acquisition and control 
scheme was used to take these data. Figure 3 is a 
flow chan showing the sequence of events for a scan. 
Two computers were used. The first computer 
controlled the movements of the arcjet and triggered 
the second computer to acquire data once the number 
of data points, the physical spacing between points, 
and the duration of each exposure were entered. The 
first computer triggered data acquisition via a TTL 
pulse. This pulse was sent to a pulse generator 
which triggered the diode array through the diode array 
controller which was interfaced with the second 
computer. The data were then saved on this computer 
before the first computer moved the thruster to the 


next location and again triggered the diode array. This 
process continued until the full scan was completed. 
After a particular scan was completed the electrical 
configuration was changed and the process was 
repeated 



Figure 3: Data Acquisition System 
Flow Chart 


Results and Discussion 
Spectra. The entire visible spectrum, for 
configuration 1, measured looking through the plume 
at the nozzle exit plane on centerline, is shown in 
Figure 4. The most prominent emissions were 
identified based on tabulated data. 15 * 16 The resolution 
of the spectroscopic system was sufficient to resolve 
the molecular band structure. The H a and Hp lines are 
not shown to scale in order to allow sufficient detail 
of the lower intensity emissions. The presence of 
atomic hydrogen, molecular nitrogen ions, and NH at 
the exit plane indicates that dissociation, ionization, 
and recombination occur in the anode/nozzle. There 
was no significant difference between this spectrum 
and spectra measured in the plume of a standard 
laboratory model arcjet with a solid tungsten anode 
operating under similar conditions. This indicates 
that the physical processes occurring within the 
plume in both cases were very similar and that the 
segmented anode arcjet accurately simulates a standard 
arcjet. However, differences may exist between this 
spectrum and the one obtained from an arcjet run on 
the actual decomposition products of hydrazine. The 
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continuum radiation between 5000 and 6500 A was 
partially the result of the dark current signal from the 
photomultiplier tube and ambient light 

Radial Intensity Profiles. The radial intensity 
variation from the Hp emission line at 4861 A for 
four different axial positions is shown in Figure 5. 
These data represent total intensities integrated over 
the entire line. At any given spatial location the 
cross sectional dimensions of the portion of the 
plume imaged on the entrance slit of the spectrometer 
was 1 micron wide by 1 mm high. Because the 
arcjet was moved 0.05 mm vertically between data 
points, these total intensities represent moving 
averages. The curves Fit to each of the radial profiles 
were calculated assuming a radial intensity 
distribution of the form I = I mix a br2 . 



Figure 5: Radial Intensity Profiles of Hjj emission at 
the Nozzle Exit Plane 


The radial intensity profiles were inverted using 
the inverse Abel transformation discussed in the 
Appendix. Figure 6 is a sample of an inverted profile 
for the H a emission line. There is a significant 
amount of scatter in the inverted values near 
centerline. Emissions from the central portion of the 
plume constituted only a small fraction of the total 
light observed along a line of sight. As a result, 
small fluctuations in the measured total intensity near 
centerline propagated through the inversion, resulting 
in larger fluctuations in the inverted intensities near 
centerline. In every case, the curve fit to the radial 
intensity profile was also inverted. Each point of the 
inverted radial intensity profile is proportional to the 
population density of that particular electronic state. 
The radial decrease in inverted intensity shown in 
Figure 6 shows how the population of this particular 
electronic state varies with radial distance. The 
populations were essentially undectactable at radii 
greater than about 4 mm. 



Figure 6: Inverted Intensity from Ha, Configuration 
2, 2.5 mm from the exit plane 


Axial Intensity Profiles. The axial variation in 
the intensity was measured along nozzle centerline for 
H a , Hp, and NH for each of the three electrical 
configurations. Figure 7 shows the axial variation 
of H a , Hp and NH A 3 n -> X 3 £ (0,0) (3361 A) for 
configuration 1. The plot has a semi-log scale and 
the observed linearity indicates the decays were 
exponential. The intensities for H a and Hp are total 
intensities. The intensities for NH(0,0) are peak 
intensities and the units are arbitrary. This is because 
themolecular spectra from the NH A 3 n — > X 3 L (1,1) 
transition at 3370 A and the bands from N 2 C 3 n 
B 3 n (0,0) and (1,1) transition at 3371.3 A and 3339 
A overlapped the NH(0,0) rotational-vibrational band, 
greatly complicating a determination of the total 
energy of any of these bands. 



Figure 7: Axial Intensity Variation of H a , Hp and 
NH(0,0) for Configuration 1 


For a particular electronic transition, the axial 
decay rate of the emitted light intensity can be 
estimated based on the rate of expansion and the 
relative rates at which the upper excited state was 
populated and depopulated. The populating 
mechanisms are electronic excitation from lower 
energy states and radiaLive decay from higher energy 
states. An estimation of the rate with which the state 
is populated due to collisional excitation required an 
estimate of the number densities of the collisional 
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partners and appropriate collisional cross sections. 
An estimation of the rate with which the state is 
populated due to radiative decay from higher energy 
states requires an estimate of the number densities for 
each of these and the appropriate transition 
probabilities. Depopulation occurs because of 
radiative decay, collisional de-excitation, and 
collisional excitation to higher energy states. 

Based on the radiative lifetimes for the transitions 
measured, the times needed for the observed intensity 
decreases were calculated assuming that the decrease in 
intensity was solely the result of radiative decay, i.e. 
no collisional de-excitation and no significant 
repopulation in the plume. The calculated times for 
the observed rate of decay for NH indicates a velocity 
approximately three times lower than for atomic 
hydrogen, if the NH molecules were expanded at the 
same rate or a slower rate than the atomic hydrogen. 
Because velocity is inversely proportional to the 
square of particle mass, higher velocities were 
expected for the lighter species. However, the times 
calculated for atomic hydrogen using the radiative 
lifetimes of Ha and Hp were significantly different. 

If the decay would soley be governed by the 
transition probability, the shorter time calculated for 
the axial decay of H a would indicate that its upper 
energy state, the third energy level, was repopulated at 
a higher rate than the upper energy state of H p , the 
fourth energy level. The rale of collisional excitation 
for these two states was estimated from the excitation 
cross sections as presented by Drawin. 17 The 
estimated rate of repopulation due to collisional 
excitation from the ground stale indicated that, while 
collisional repopulation of the third energy level 
occurred more frequently than repopulation of the 
fourth energy level, significant collisional 
repopulation of these states was unlikely. Therefore, 
repopulation due to radiative decay of higher energy 
states was the significant repopulating mechanism. 

Based on a consideration of the transition 
probabilities and oscillator strengths for various other 
hydrogen transitions 15 it was likely that third and 
fourth energy levels were populated at comparable 
rates from decay of the fifth energy level and above. 
However, the preferential repopulation of the third 
energy level could occur as a result of the P a 
transition, from energy level 4 to 3. This transition, 
which is not visible (18751 A), accounted for the 
observed rate of decay for H a . Since the transition 
probability for P a was smaller than that for H a , the 
dominant rate controlling step in the decay of energy 
level 3 was the rate at which it was repopulated from 
the fourth energy level. Calculations of the time 
assosciated with the intensity decrease of H a using the 
decay rate of P a agreed with the values calculated for 


Hp, because the radiative lifetime for P a is close to 
that of Hp. 

Electron Number, Density and Electronic 
Excitation Temperature. The determination of the 
radial distribution of the electron number density from 
the full width at half maximum of the Hp line due to 
Stark broadening was not possible using the data 
taken with the linear diode array due to instrumental 
broadening. For all cases the measured FWHM was 
0.4 A indicating narrower lines could not be resolved. 
Calculations verified that this was the result of the 
resolution of the linear diode array. Although the 
limited resolution of the system precluded a 
determination of local electron number densities, a 
FWHM of 0.4 A permitted a calculation of an 
absolute upper limit for the electron number density. 
As discussed previously, to obtain the electron 
number density from a line profile, the line must be 
deconvoluted to determine the amount of Doppler and 
Stark broadening. To determine the amount the line 
was broadened due to Doppler broadening, a kinetic 
temperature for electronically excited hydrogen of 
5000 K, based on measurements to be discussed in 
the following section, was used. Using this 
temperature, the calculated electron number density 
for an Hp line whose FWHM is 0.4 A was 2.2xl0 13 
cm' 3 . Thus the electron number densities at each 
point measured were below this value. 

Subsequent to this test, the FWHM for the Hp 
line was measured in the plume of an identical 
laboratory-model thruster with a solid anode operating 
at the same conditions. These later measurements 
were taken using the photomultiplier tube. As a 
result of the higher resolution of this detector, the 
instrumental broadening was reduced and the Hp line 
was resolved. The FWHM measured at the exit 
plane on centerline was 0.34 A. These 
measurements were non-inverted, line-of-sight 
intensities values. For a temperature of 5000 K, this 
FWHM corresponds to an electron number density of 
1.8xl0 13 cm 3 . 

As noted in the Analysis section, if an 
equilibrium or Boltzmann distribution exists between 
the different electronically excited states of an atom, 
the line intensity method can be used to determine the 
electronic excitation temperature. Temperature 
profiles were calculated using the inverted intensities 
from H a and Hp emission lines assuming 
equilibrium. For electrical configurations 1 and 2 
flat temperature profiles were obtained with 
temperatures of approximately 5000 K and 5500 K. 
While these temperatures seem reasonable, the 
intensities measured for electrical configuration 3 
produced negative values for electronic excitation 
temperature. 
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Based on calculations of the rate of coliisional de- 
excitation from the excitation cross sections, 17 
radiative decay was the dominate de-excitation 
process. Because the radiative lifetime of H a is 
approximately five times less than that for Hp, one 
would expect the upper excited state of Hq should be 
depopulated more rapidly than Hp. For this case, the 
line ratio method would yield falsely high 
temperatures for a slight deviation from equilibrium 
and negative values for more substantial deviations. 
The ratio of the intensities from the two transitions is 
shown in Figure 8. The increasing intensity ratio 
with distance indicated a more rapid depopulation of 
Ha- 



Figure 8: Axial Hp/H a Intensity Ratio 

For an assumed electron temperature of 5600 K 
(0.5 eV), electron number densities greater than 10 14 
cm* 3 are required to satisfy the conditions for PLTE 
presented in the analysis section. This is 
approximately an order of magnitude higher than the 
electron number density estimated from the observed 
Hp line width. The negative temperature values 
calculated from the H a and Hp line intensities for 
configuration 3 and the measured electron number 
density indicate a non-Boltzmann energy distribution 
between these upper excited states. As a result, the 
line intensity ratio is not valid for the determination 
of excitation temperatures for any of the 
configurations. 


was estimated to be ±50 K. No conclusive 
evaluation of the rotational temperature from the NH 
A 3 n — > X 3 I (0,0) and (1,1) transitions was possible 
because of the overlap of these two bands and the two 
N 2 bands in this region mentioned previously. 
However, a temperature near 2500 K , measured for 
NH from the C l U -> B 1 !, transition was likely. 

Table IV. Electronic Transitions and Rotational 


Temperatures 


Electronic Transition 

X, A 

Troll K 

n 2 c 3 n„-> B 3 n g (0,2) 

3804.9 

800 

n 2 c 3 rr u -> B 3 n g (i,3) 

3755.4 

700 

N 2 + B 2 ^-* X%(0,0) 

3914.4 

1800 

N 2 * B 2 2^— >X 2 Zg((l,l) 

3884.3 

1700 

NH C'n -> B 1 ! (0,0) 

4504.0* 

2500 

NH A 3 n -> X 3 L (0,0) 

3360 

— 

nh A’n-^^in,!) 

3371 

_ 

* Q head, R head at 4525 A and P head at 4538 A 


Because the relaxation times for rotational 
equilibration are of the same order as that for 
adjustment of the translational motions of 
molecules, 18 * 19 (i.e. only a few collision are needed), 
in the absence of local thermodynamic equilibrium, 
the rotational temperature may provide a good 
estimate of the kinetic temperature of a molecule. 
However, if coliisional rates are insufficient to 
provide equilibrating rotational energy transfer, the 
molecules are rotationally frozen and the rotational 
temperature reflects the kinetic temperature at the 
point the molecules were frozen, No attempt was 
made to estimate the collision frequencies for the 
molecular constituents due to the difficulty in finding 
accurate estimates of the coliisional cross sections, 
however, the relatively low rotational temperature 
values determined for N 2 at the nozzle exit plane 
compared to estimates of the rotational temperature in 
the constrictor suggest the nitrogen was not 
rotationally frozen. Therefore, the temperature of 
approximately 750 K reflects the degree to which the 
nitrogen was cooled in the gasdynamic expansion. 


Rotational Temperature. Molecular spectra from 
N 2 , N 2 +, and NH were used to determine rotational 
temperatures. The thruster was operated in 
configuration 1 for all of these data. The rotational 
temperatures determined from the listed electronic 
transitions and the location of the band heads are 
listed in Table IV. The intensities of these molecular 
spectra were measured with the photomultiplier tube. 
These values were not inverted due to the added 
complexity of inverting an entire molecular band 
prior to the temperature analysis. Therefore, the 
intensities are line of sight averages across the nozzle 
exit plane on centerline. As a result, the temperatures 
are mean values for the plasma in this region. Based 
on the radial intensity profiles, the signal was 
averaged over approximately 5 mm. The uncertainty 


The rotational temperature for the molecular 
nitrogen ions was determined to be approximately 
1000 K hotter than the molecular nitrogen. Because 
the electron - ion coliisional cross sections are much 
larger than the electron - neutral coliisional cross 
sections, heating via electronic collisions occurred 
more frequently for N 2 + . Therefore, the rotational 
energy distribution for N 2 + reflected a higher 
temperature than for N 2 . The rotational temperature 
for the NH molecules was substantially higher than 
the values for N 2 and N 2 + . As previously indicated, 
the NH molecules resulted from recombination of 
dissociated hydrogen and nitrogen within the nozzle 
This recombination could occur as the result of two 
body collisions or three body collisions involving 
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electrons, heavy panicles, or the nozzle walls. The 
rates of radiative recombination, (i.e. two body 
recombination), and three body recombination were 
estimated based on recombination coefficients 
presented in Mitchner and Kruger. 20 These 
calculations indicate the dominant recombination 
mechanism was three body recombination. In this 
case, some or all of the recombination energy was 
carried off as increased kinetic energy of the 
collisional partner or as heating of the nozzle walls. 
The remainder of the recombination energy remained 
with the recombined molecule. With sufficient 
vibrational-rotational energy transfer, a portion of the 
recombination energy increased the molecule's, 
rotational temperature. Therefore, the recombined NH 
molecule was rotationally hotter than the other 
molecular constituents. 

Electrical Configuration. The electrical 
configuration had a measurable effect on the overall 
intensities of the observed H a , Hp, and NH 
transitions. For example, the inverted H a radial 
intensity profiles from the nozzle exit plane for each 
configuration are shown in Figure 9. The intensities 
at configuration 2 were more than three times brighter 
than those measured for configuration 1 and monT 
than 20 times brighter than those for configuration 3, 
Similar results were obtained at other axial positions 
for the H a , Hp, and NH A 3 n -» X 3 E transitions. The 
lower intensities for configuration 1 and configuration 
3 are the result of a lower population of the excited 
states. This suggests that significant population of 
excited states is occurring only within the current 
attachment regions of the anode. 



Figure 9: Inverted H a Intensities at the Exit Plane for 
Each Electrical Configuration 


The changes in the measured intensity for the 
different electrical configurations were attributed to 
changes in the rate of electron collisional excitation 
within the nozzle, which determines the population of 
excited states. The rate of collisional excitation 
depends on the number density of the species being 
excited, the electron number density, the electron 
energy distribution, and the appropriate collisional 
cross sections. There were significant changes in the 
distribution of the current carrying electrons in the 
nozzle with configuration. A summary of the current 
distribution and overall power for each configuration 
is shown in Table IV. Based on spectroscopic 
measurements inside the anode of a similar 
laboratory-model arcjet the electron number density 
inside the nozzle decreases from approximately 10 15 
cm -3 near the constrictor to approximately 10 14 cm* 3 
near the nozzle exit, 21 However, the number of 1 eV 
electrons needed to carry 10 Amperes of current is 
approximately 10 u cm’ 3 . Thus, changes in the 
distribution of current carrying electrons likely had a 
small effect on the overall electron number densities 
inside the nozzle. 

Since the electron number density does not 
change significantly with electrical configuration, it 
is reasonable to assume that changes in the electron 
energy distribution function were responsible for the 
observed changes in the excitation spectrum. For 
example, by changing from configuration 2 to 
configuration 3, the electron acceleration due to the 
electric field is forced to occur in a higher pressure 
region of the nozzle. Thus, more collisional energy 
exchange between the electrons and heavy particles 
occurred in configuration 3 limiting the number of 
energetic electrons available for collisional excitation. 
In configuration 2 the electron acceleration due to the 
electric field continues downstream to segment 5. In 
this region there was less collisional energy exchange 
with heavy particles due to lower total number 
densities. Consequently, the electrons were 
accelerated to higher velocities than in configurations 
1 and 2 and, more energetic electrons were available 
for collisional excitation 

In configuration 3, the electrons were no longer 
accelerated due to the electric field in the region of the 
anode corresponding to segment 5. Therefore, in this 
region the electron collisional excitation depleted the 
high energy electrons. Electron self equilibration 


Table IV: Voltage/Current Characteristics 


Segment 1 

Segment 2 

Segment 3 

Segment 4 

Segment 5 

voltage 

current 

voltage current 

voltage current 

voltage 

current 

voltage 

current 

Configuration 1 54 

0.0 

135 2.0 

135 1.9 

135 

1.9 

135 

4.3 

Configuration 2 55 

0.0 

100 0.0 

no o.o 

125 

0.0 

140 

10.0 

Configuration 3 138 

0.4 

138 4.4 

138 4.4 

138 

0.8 

122 

0.0 



Voltages in Volts; Currents in Amperes 
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rapidly seeks to reestablish a Maxwellian energy 
distribution, however, at a slightly reduced 
temperature. A lower electron energy distribution 
further reduces the rate of electron collisional 
excitation. It seems that this electron cooling 
mechanism occurred with such rapidity that the rate of 
collisional excitation quickly approached the rate of 
radiative decay within this part of the nozzle. This 
resulted in the more pronounced deviation from 
equilibrium for the excited states of atomic hydrogen 
indicated by negative temperature values calculated 
using the line intensity ratio for configuration 3. 

Frozen How Losses. In the arcjet, electrons are 
accelerated by the electric field and imparl energy to 
the propellant via collisional processes. Some of this 
energy is invested in excitation, ionization and 
molecular dissociation. The energy that remains in 
these internal modes at the nozzle exit is called the 
frozen flow loss. Based on the overall device 
efficiency and reasonable estimates of other loss 
mechanisms it is likely that frozen flow losses 
account for 30 to 40 percent of the total energy input 
to the device. 22 

The variation in the measured intensities 
demonstrated that the location of the current 
attachment on the anode significantly effects excited 
state number densities in the plume. Results reported 
in a companion paper indicate that the location of the 
current attachment on the anode had a negligible effect 
on the measured thrust of the device. 10 These results 
imply that the frozen flow losses due to electronic 
excitation were negligible. Because 12. l.and 12.7 eV 
are required to excite the third and fourth energy levels 
of hydrogen, the upper excited states of H a and Hp, 
this result has significant implications about other 
frozen flow losses in kilowatt-class arcjets. 

Consider the ionization and dissociation energies 
shown in Table V. Recent measurements indicate 
that the electron temperature in the anode is below 2 
eV. 21 If the electron energy distribution for this 
temperature is Maxwellian, in the absence of stepwise 
excitation, the electrons in the high energy tail were 
responsible for ionization and the excitation of upper 
excited states with energies approaching the 
ionization potential, such as H a and Hp. The fraction 
of electrons in this high energy tail is small, (eg. the 
fraction of electrons having energies above 13.6 eV is 
0.005 for a 2 eV electron temperature). Given that 
significant changes in the frozen flow losses due to 
excitation of hydrogen were insignificant from a 
performance standpoint, it is unlikely that other 
processes requiring these energetic electrons were 
dominant. Therefore, ionization, although integral 
to the operation of the arcjet, was responsible for 
only a small part of the frozen flow losses, because of 
the high energy thresholds relative to the electron 


temperature. Ionization losses of approximately 10 
Watts, or less than 1 percent of the total energy 
input, were estimated. This estimate was based on 
the assumption of quasi-neutrality, (ie. the ion 
number density approximately equal to the electron 
number density), an ion exit velocity of 5000 m/s, 
based on the measured specific impulse, and an 
average ionization energy of 14 eV. For ionization 
losses on the order of 10 Watts, less than one percent 
of the flow was ionized. This agrees with the 
ionization fraction reported by Zana. 5 


Table V. Ionization and Dissociation Energies 


Species 

Ionization *, eV 

Dissociation , eV 

H 

13.60 

— 

h 2 

15.48 

4.53 

N 

14.53 

— 

n 2 

15.58 

9.85 

NH 

13.10 

3.26 

___ 

* first ionization state 


Therefore, because frozen flow losses due to 
electronic excitation and ionization were small, the 
dominant frozen flow loss mechanism was 
dissociation. This is because the energy thresholds 
for dissociation are much closer to the electron 
temperature. As a result, the number of electrons 
capable of collisional dissociation is relatively large, 
(e.g. the fraction of electrons with energies above the 
4.5 eV threshold for dissociation of hydrogen is 0.2 
for a 2 eV electron temperature). A smaller number 
of electrons have energies above the 9.9 eV 
dissociation energy of nitrogen. Minimal amounts of 
atomic nitrogen were detected at the exit plane. The 
presence of excited NH in the emission spectra 
indicated that significant amounts of atomic nitrogen 
had recombined in the nozzle. However, due to the 
small dissociation energy of the NH molecule, this 
recombination recovered only some of the energy lost 
due to dissociation. A total of over 14 eV is required 
to dissociate both a hydrogen and a nitrogen 
molecule; however, only 6.5 eV are recovered by two 
NH recombinations. 

Conclusions 

Light intensities resulting from the de-excitation 
of electronically excited species in the near field of the 
plume were measured. Excitation temperature, 
rotational temperature, and electron number density 
were evaluated based on these intensity 
measurements. The effect the anode current 
distribution had on the characteristics of the plume 
was investigated. 

The complete emission spectrum, measured from 
3200 A to 7200 A at the nozzle exit plane detected the 
electronically excited species N 2 , N 2 + , NH, and H, 
indicating excitation, dissociation, ionization and 
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recombination in the nozzle. Radial intensity 
profiles, spaced 0.05 mm apart, were measured for the 
H a , Hp, and NH A 3 n X 3 L (0,0) transitions at four 

different axial locations. These line of sight 
intensities measurements were deconvolved to give 
the radial intensity distribution using an inverse Abel 
transformation. Axial intensity profiles taken on 
centerline indicated an exponential decay in excited 
state population for Hq, Hp, and NH. The rate of 
axial decay indicated lower velocities for NH than H 
in the plume and population of the third excited 
energy state of hydrogen from the decay of higher 
energy levels. The suitability of a line intensity ratio 
method, which requires a Boltzmann energy 
distribution between excited states for the 
determination of an excitation temperature, was 
evaluated. In some instances, this method resulted in 
negative temperatures. This indicated a non- 
equilibrium energy distribution between excited 
energy states of hydrogen in the plume. This was 
analytically verified. Rotational temperatures 
evaluated based on the emission spectra for several 
molecular transitions ranged from 750 K for N 2 to 
2500 K for NH, Based on these results, the arcjet 
plume was found to be a highly non-equilibrium 
plasma. 

Anode electrical configuration was found to have 
a large effect on the spectral intensities measured in 
the plume. These data indicate that significant 
population of excited states occurs mainly in the 
current attachment regions in the nozzle. 
Therefore, by forcing the arc to attach farther 
upstream, frozen flow losses losses due to electronic 
excitation were significantly reduced. The reduction 
in the frozen flow losses due to excitation was 
determined to be negligible from a performance 
standpoint based on thrust measurements reported in a 
companion paper. 10 As a result, the dominant frozen 
flow loss mechanism in kilowatt class arcjets run on 
hydrogen/nitrogen mixtures was dissociation. 
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Appendix 

Abel inversion The spectral intensities measured 
in this investigation were line of sight values, 
measured by looking through the plume from the 
side, perpendicular to the arcjet center line. For line 
of sight measurements, the intensity values 
correspond to integrated values of intensity as a 
function of position. If the plume is axisymmetric, 
the emission coefficient at any point, £, is only a 
function of the radial coordinate, r. This situation is 
depicted schematically in Figure 10. 



Figure 10: Projection of an Axisymmetric Function 

If there is no significant reabsorption of the 
emitted radiant energy in the plume, it is optically 
thin. For an axisymmetric, optically-thin plume the 
observed intensity at a location x is given by: 

I(x) = J e(r) dy = J e(VA?) dy (Al) 

There are several methods to determine the 
emission coefficient, e(r), from Equation (Al). The 
methods for reconstructing this radially dependent 
emission coefficient from the projected intensity 
values are known as inverse Abel transformations. 
Several techniques have been used to perform inverse 
Abel inversions. The method presented here is based 
on that originally presented by Sudharsanan. 23 The 
algorithms utilized and the peak finding routine were 
independently derived and are presented here. The 
Fourier transform of Equation (Al) results in the 
following expression: 

/ °° -j2Kfx 

I(x)e dx (A2) 

By changing to polar coordinates and substituting 
for the zeroth order Bessel function, it can be shown 


that equation (A2) is equal to the zeroth order Hankel 
transform of e(r). 23 Therefore, the emission 
coefficient can be obtained from the inverse Hankel 
transform. 



G(0fJo(27rtr)cf 


(A3) 


In actual practice the intensities are only 
available at discrete positions and the center line of 
the data is not known. This is accommodated using a 
discrete Fourier representation of the projected 
intensity values. The Fourier transform can be 
evaluated using the following direct Fourier transform 
(DFT) algorithm: 

Gi=P(')exp(— ) (A4) 

i ■ 0 


In order to evaluate the emission coefficient, the 
inverse Hankel transform is performed. Since the 
intensity profile is axisymmetric, the inverse Hankel 
transform need only be calculated for half the 
intensity profile. However, to do this the transformed 
intensity data must be shifted such that the axis of 
symmetry is at the origin. For noise-free data 
correctly centered on the interval i = 0 to N-l, the 
imaginary part of the DFT is zero. The magnitude of 
the shift needed to properly center the data can be 
determined by minimizing die imaginary component 
of the DFT of the actual data. The values calculated 
using the algorithm in equation (A4) are then shifted 
the appropriate amount using the DFT shift 
theorem 24 to place the center line at i = 0. The 
emission coefficient is then evaluated using the 
following discrete representation: 

«®-r 1L ilWy) (A5) 

N (Ax) l-o 

Figure 11 shows the inverted results for two 
analytic test functions. Figure 11(a) shows inverted 
values for an ellipse. The exact solution is a step 
function. There is ringing in the calculated solution 
at the singularity. However, this type of profile 
would not likely be encountered and is included only 
as a worst case. The solution shown in Figure 1 1(b) 
is for a Gaussian input function. The exact solution 
and the calculated values are coincident. This is 
closer to the cases expected in emission spectroscopy 
experiments. 
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Figure 11: Verification of Inverse Abel 
Transformation with Analytic Functions 
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electron number density of 2x 10 13 cm -3 at the exit plane was determined based on Starlc broadening of the 
line. Rotational temperatures of 750 K, 1750 K, and 2500 K were determined for N 2 , N 2 * , and NH respectively. 
The results of this investigation demonstrate that the location of the current attachment on the anode has a 
measurable effect on the electronically excited species in the plume and that dissociation is the dominant frozen 
flow loss mechanism in low power arcjets. 


- 

17. Key Words (Suggested by Author(s)) 
Arcjet 

Spectroscopy 

Plume 

18. Distribution Statement 

Unclassified - Unlimited 
Subject Category 20 

19. Security Classif. (of this report) 
Unclassified 

20. Security Classif. (of this page) 
Unclassified 

21. No. of pages 

17 

22. Price* 

A03 


NASA FORM 1626 OCT 86 


*For sale by the National Technical Information Service, Springfield, Virginia 22161 




